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SYMBOLS 


A — Azimuth, measured from the direction of mo- 
tion of the airplane (in aircraft fire) 

A — Altitude function (Eq. 63) 

A — Axial moment of inertia 
B — Drag coefficient in lb/in 2 -ft (Eq. 19) 

C — Ballistic coefficient (Eq. 22) 

D — Drag 

D — Angular deflection 
F — Temperature in degrees Fahrenheit 
G — Drag function (Eq. 20) 

H — Ratio of air density at altitude y to that at 
the surface (Eq. 48) 

H — Height of target 
I — Inclination function (Eq. 62) 

K — Coefficient (See Eq. 16 and 34) 

L — Cross wind force (lift) 

M — Mach number (Eq. 18) 

M — Pseudo-Mach number (Eq. 65) 

M — Overturning moment about the center of 
gravity 

M — log 10 e=0.43429 (approx.) 

N — Spin 

P — Performance parameter (Eq. 27) 

P — Siacci range : the distance measured along the 
line of departure to a point directly above 
the projectile (Eq. 54) 

Q — Drop of the projectile (Eq. 55) 

R — Reynolds number (Eq. 17) 

R — Range 

S — Space function (Eq. 60) 

T— Total time of flight 
T- — Time function (Eq. 61) 

U— Upper limit of integration of Siacci functions 
V — Velocity 
W — Wind speed 
X — Total range of a trajectory 
Y — Altitude of an airplane above sea level 


Z — Linear deflection 
Z — Zenith angle 
i — Velocity of sound 

Diameter of the largest cross section of the 
projectile (in Gavre drag law, Eq. 23) 
c' — Damping coefficient due to yawing moment 
and crosswind force (Eq. 51) 
c" — Damping coefficient due to drag (Eq. 33) 
d— Caliber, or maximum diameter of projectile 

f — Drag coefficient (in the Gavre drag law, Eq. 
23) 

g — Gravitational acceleration 
h — Logarithmic rate of decrease of air density 
with altitude : 0.000,031,58 per ft 
h — Yawing moment damping factor 
i — Form factor (Eq. 21) 

k — Logarithmic rate of decrease of sonic velocity 
with altitude 0.000,003,01 per ft 

k — Ratio of American to French drag function 
(Eq. 24) 

m — Mass of the projectile 
m — Slope of the tangent to the trajectory 
n — Logarithmic rate of change of muzzle velocity 
with projectile weight 
p — Space function (Eq. 28) 
p — Projectile weight 
q x — Inclination function (Eq. 30) 
q — Altitude function (Eq. 31) 
s — Stability factor 

t — Time of flight (sometimes the symbol t f is used 
to distinguish the time of flight from the 
Siacci time function) 

t — Time function (Eq. 29) ^ " 

u — Velocity of the projectile relative to the air 

u — Component along the line of departure of the 
velocity relative to the air : the argument of 
the Siacci functions 

v — Velocity of the projectile relative to an inertial 
system 


IV 



w — True air speed of the airplane in aircraft fire 
x — Horizontal distance from the origin of the 
trajectory 

y — Vertical distance from the origin of the tra- 
jectory 
z — Drift 

A (capital delta) or A' — Increment of (followed 
by a variable) 

A — Density of air (in Gavre drag law, Eq. 23) 

Y (gamma) — Factor inversely proportional to the 
ballistic coefficient (Eq. 66) 

8 (delta) — Yaw 

e (epsilon) — Angle of site (Eq. 72) 

$ (theta) — Angle of inclination of the trajectory 
x (kappa) — Crosswind force damping factor 
(mu) — Viscosity of the air 
p (rho) — Density of the air 
t (tau) — Temperature of the air 
<p (phi) — Angle of departure 
G> (omega) — Angle of fall 
Superscripts 

One dot (example: x) — Derivative with respect 
to time. 

Two dots (example: y) — Second derivative with 
respect to time. 

Subscripts 

0 — At t = 0, y = 0, or 8 = 0. 


1 — Pertaining to, or relative to, Projectile Type 1 

2 — Pertaining to, or relative to, Projectile Type 2 
8 — Pertaining to, or relative to. Projectile Type 8 
.1— First revision 

.2 — Second revision 

C — Due to 1 % increase in ballistic coefficient 
D — Drag 

H — Due to height of target 
R — Range component of 
e — Effective 

p — Due to projectile weight 
s — Summital : at the summit of a trajectory, where 
y = 0 

s— Standard : pertaining to standard conditions 
t — Pertaining to, or relative to, some Typical Pro- 
jectile 

v — Due to 1 unit increase in muzzle velocity 
w — Due to wind 
z — Cross component of 
8 — Yaw 

x — Due to increase in the temperature of the 
atmosphere 

y — Due to an increase in angle of departure 

( 0 — Terminal : where y = 0 on the descending 
branch of the trajectory 
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PREFACE 

This pamphlet, one of the series which will comprise the Ordnance Engineering 
Design Handbook, deals with the subject of Exterior Ballistics, covering Tra- 
jectories, Differential Effects and Data for Projectiles. The pamphlet pertains 
to that phase of projectile life starting at the gun muzzle (or, for a rocket, at 
burnout) and ending at the point of burst or impact. The subject is covered in a 
manner intended primarily to assist the ammunition designer in recognizing and 
dealing with the design parameters of projectile weight, shape, muzzle velocity, yaw, 
drag and stability, in designing for optimum results in range, time of flight and 
accuracy. It is hoped that it will also prove helpful to other personnel of the Ord- 
nance Corps and of contractors to the Ordnance Corps. Because of its complexity, 
it is possible to cover the subject only briefly in this volume, but an attempt has been 
made to incorporate references which will permit the designer or the student to 
explore in more detail any phase of the subject. 

The ballistics of bombs and of intermediate and long range missiles are not 
covered in this pamphlet. Neither are the problems associated with Interior 
Ballistics or of Terminal Ballistics. These subjects will be covered by other 
pamphlets. 
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TRAJECTORIES, DIFFERENTIAL EFFECTS, AND DATA 

FOR PROJECTILES 


1. Introduction to Exterior Ballistics 

Ballistics is the science of the motion and 
effects of projectiles. Exterior Ballistics deals 
with that part of the motion between the pro- 
jector and the point of burst or impact. 

The projector may be any sort of weapon: a 
sling-shot, a bow, a mortar, a rifle, a rocket 
launcher, an airplane, etc. A projectile may be 
any body that is projected by such a weapon. 
However, in this pamphlet, the term “projectile” 
will be restricted to those missiles that are in free 
flight through the air after being propelled by a 
gun or a rocket motor. This term does not in- 
clude bombs, which are released from airplanes; 
they will be considered in another pamphlet. 

The motion of a missile in flight is really 
quite complicated. Even a fin-stabilized pro- 
jectile or bomb that is not spinning about its 
longitudinal axis is swinging about a transverse 
axis, and this affects its motion. The nose of a 
spinning shell precesses about the tangent to its 
trajectory, and its attitude gives rise to forces that 
not only affect its motion in a vertical plane but 
also move it out of this plane. To simplify the 
calculations, however, it is usually assumed that 
a missile acts like a particle subject to two forces: 
(1) gravity, which is directed vertically down- 
ward, and (2) drag, which is directed opposite to 
the direction of motion. Gravity is assumed con- 
stant : until 1956, the gravitational acceleration was 
usually taken to be 32.152 ft/sec 2 (9.80 m/sec 2 ) ; 
and was used for all tables referred to in this text ; 
now, the standard value is 32.174 ft/sec 2 (9.80665 
m/sec 2 ). 1 The drag depends on the size of the 
missile, its velocity, the density of air, and other 
factors. Under these conditions, the motion fol- 
lows a smooth, planar path, which is called an 
“exterior ballistic trajectory,” or “trajectory” for 
short. 

The initial conditions for the trajectory of a 
particular missile are its initial velocity and angle 
of departure.* For a projectile, the initial velocity 

♦The term angle of departure and the term elevation 
as used in this pamphlet represent the angle from the 
horizontal to the line of departure and therefore would be 
more accurately and completely termed quadrant angle of 
departure . 


is the muzzle velocity and angle of departure is 
approximately’ the elevation of the weapon 
(fictitious values of these quantities may be used 
in the case of a rocket that burns outside the 
launcher). A bomb is released with the speed of 
the airplane, generally in a horizontal direction. 
If a gun is mounted in an airplane or other 
moving vehicle, the initial conditions are affected 
by the motion of the vehicle. In any case, it is 
customary to compute trajectories for a few se- 
lected values of initial velocity and elevation. 

Certain standard conditions are specified for a 
normal trajectory : not only the initial velocity and 
angle of departure, but also the weight of the 
projectile and a factor that depends on its shape ; 
the density, temperature, and lack of motion of 
the air, and the height of the target relative to the 
projector. In order to determine the effects of 
variations from these conditions, it may be con- 
venient in some cases to compute non-standard 
trajectories; however, they can usually be esti- 
mated in other ways, as will be explained later. 

2. Vacuum Trajectories 
a. Equations of Motion 

A rough idea of an exterior ballistic trajectory 
may be gained by neglecting the air resistance, 
considering only the effects of gravity. Then if 

x is the horizontal distance from the origin, 
y is the vertical distance from the origin, . 
t is the time of flight, 
g is the gravitational acceleration, 

and dots denote derivatives with respect to time, 
the differential equations of motion are 

X = 0, y = — g. (1) 

The initial conditions, denoted by the subscript 0, 
are 

x» = v„ cos •a,, y„ — v„ sin a o , x« = 0, y« = 0, 

where v is the velocity, 
a the angle of inclination. 

By integrating equations (1), we find that, at 
any time, 

x = v 0 cos a„, y = v 0 sin a, — gt, (2) 

and by integrating (2), 


1 
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x = (v 0 cos *o)t, y = (v 0 sin $ Q )t — gt 2 /2. (3) 
b. Characteristics of Trajectories 

From Equations (3), we find that the time at 
any horizontal range is 

t^x/VoCOS'&o (4) 

and hence the ordinate is 

y = x tan # 0 — gx 2 / 2v 0 2 cos 2 #<,• (5) 

This is a parabola with a vertical axis. 

The total range X is found by letting y —0 in 
(5) and solving for x: 

X = (2v 0 2 /g) sin *8 0 cos £ 0 = (v 0 2 /g) sin 2'8 0 . (6) 
This shows that the range is a maximum when 
$o = 45°. 

The total time of flight T is found by substi- 
tuting (6) for x in (4) : 

T = (2v 0 /g) sin^ 0 . (7) 

The components of the terminal velocity are 
found by substituting (7) for t in (2) : 

X = Vo cos 8 0 , Y = — v 0 sin £ 0 . (8) 

Hence the terminal velocity is 
v<o = v 0 , (9) 

and the angle of fall (considered positive) is 

<*> = 8 0 . ( 10 ) 

By differentiating equation (5) with respect to 

x, we find that the slope of the trajectory at any 
point is 

~~ = tan — gx/v o 2 cos 2 8 0 . (11) 

At the summit, the slope is 0 ; hence the horizontal 
range to the summit is 

x B = (voVg) sin # 0 cos #o = X/2. (12) 

By substituting this for x in (5), we find the max- 
imum ordinate 

y 8 = (v 0 2 /2g) sin 2 $ 0 — (1 — cos 2£ 0 ) v c 2 /4g 
= (g/8) T 2 . (13) 

By substituting (12) for x in (4), we find the 
time to the summit 

t 8 = (v 0 /g) sin * 0 = T/2. (14) 

3. Air Resistance 

a. Forces Acting on a Missile Moving in Air 
A missile moving in air encounters other forces 
besides gravity. The principal components of the 
force produced by motion through air are drag, 
crosswind force, and, for a spinning projectile, 
Magnus force. 

( 1 ) Drag is the component in the direction op- 


posite to that of the motion of the center of 
gravity. This is the only component acting on a 
non-spinning projectile that is trailing perfectly. 

(2) Crosswind Force in the component per- 
pendicular to the direction of motion of the center 
of gravity in the plane of yaw. The yaw is the 
angle between the direction of motion of the 
center of gravity and the longitudinal axis of the 
missile. This component is the cause of the drift 
of a spinning shell. 

(3) Magnus Force acts in a direction per- 
pendicular to the plane of yaw. This is the force 
that makes a spinning golf ball or baseball swerve. 

b. Simplified Assumptions 

The motion caused by these forces is explained 
in detail in the Ordnance Corps Pamphlet, “De- 
sign for Control of Flight Characteristics.” 2 In 
general, the motion of the center of gravity is 
affected by the motion about the center of gravity 
and occurs in three dimensions. To simplify the 
calculations, it is usually assumed that the only 
forces acting on the missile in free flight are 
drag and gravity. The crosswind force is implied 
in the drift, but this is determined experimentally 
as a variation from the plane trajectory. This 
approximation is satisfactory except at extremely 
high elevations, where the summital yaw is very 
large. 

c. Drag Coefficient 

The drag encountered by a projectile moving in 
air is defined by the equation of motion 

dv 

m-£-— — D — mg sin £, (15) 

where 

m — is the mass of the projectile, 
v — the velocity of the projectile relative to an 
inertial system, 
t — the time, 

D — the drag, 

g — the gravitational acceleration, 

*8 — the angle of inclination. 

The drag coefficient is defined by the relation 
K d = D/pd 2 u 2 , (16) 

where 

K d — is the drag coefficient, 
p — the air density, 
d — the caliber or maximum diameter, 
u — the velocity of the projectile relative to 
the air. 

The drag coefficient is a function of the Reynolds 


2 
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sign for Control of Flight Characteristics.” 2 In 
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calculations, it is usually assumed that the only 
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drag and gravity. The crosswind force is implied 
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c. Drag Coefficient 

The drag encountered by a projectile moving in 
air is defined by the equation of motion 

dv 

m-£-— — D — mg sin *S, (15) 

where 

m — is the mass of the projectile, 
v — the velocity of the projectile relative to an 
inertial system, 
t — the time, 

D — the drag, 

g — the gravitational acceleration, 

*8 — the angle of inclination. 

The drag coefficient is defined by the relation 
K d = D/pd 2 u 2 , (16) 

where 

K d — is the drag coefficient, 
p — the air density, 
d — the caliber or maximum diameter, 
u — the velocity of the projectile relative to 
the air. 

The drag coefficient is a function of the Reynolds 
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number R and the Mach number M. The 
Reynolds number is defined as 

R = udp/jJL (17) 

where p, is the viscosity of the air, and the Mach 
number as 

M = u/a, (18) 

where a is the velocity of sound in the surround- 
ing medium.* The viscosity is related to the 
friction of the air; the velocity of sound, to its 
elasticity. At extremely low velocities, the 
Reynolds number causes greater variation, but 
the effect of the Mach number is larger when it 
exceeds 0.5. The angle between the direction of 
motion of the center of gravity and the axis of 
symmetry of the shell, called yaw, also affects the 
drag coefficient, but this is usually stripped out of 
the resistance firing data. 

Figure 2 (p. 22) is a graph of drag coefficient 
vs Mach number for the four typical projectiles 
whose shapes are shown in Figure 1 (p. 21). 
These projectiles are described and the basic 
firings explained in Section 4. It should be noted 
that subscripts are used to distinguish the typical 
projectiles to which the drag coefficients pertain; 
a second subscript is sometimes added after a 
decimal point to denote a revision. It is evident 
that all the drag coefficients increase rapidly in 
the vicinity of Mach number 1. The main differ- 
ence between curves is in the ratio of the maxi- 
mum, which occurs between Mach numbers 1.0 
and 1.5, to the minimum subsonic value. Although 
the drag coefficient has been determined and tab- 
ulated for several other projectiles, these four 
represent the principal shapes for artillery pro- 
jectiles. 

K d is a non-dimensional coefficient, and a con- 
sistent set of units should be used in Equation 
(16). However, it is customary to express the 
density as a ratio, the caliber in inches, and the 
velocity in feet per second; then, since the drag 
is in poundals, the drag coefficient is expressed in 
pounds per square inch per foot and is denoted by 
B. Until 1956, the standard air density was 
taken as 

p 0 = 0.07513 lb/ft 8 = 5.217 X 10“ 4 lb/in 2 -ft * 

The practical drag coefficient, 

B 5.217 X 10- 4 K D , (19) 

for Projectile Types 1, 2 and 8 is tabulated as a 
function of velocity for the standard velocity of 
sound, 

a 0 = 1120.27 fps* 

* In 1956, the Department of Defense adopted the stand- 
ard atmosphere of the International Civil Aviation Organ- 
ization. 1 In this atmosphere, at the surface, the density is 


The non-dimensional drag coefficient K D for the 
90-mm HEAT Shell T108 is tabulated as a func- 
tion of Mach number. (See References 3, 4, 5 and 
6, which may be obtained from the Ballistic Re- 
search Laboratories, Aberdeen Proving Ground, 
Maryland) . 

d. Drag Function 

The drag function G is computed by the formula 
G = Bu (20) 

with u in feet per second. For convenience in 
computing trajectories, this has been tabulated 
as a function of u 2 /100 with u in meters per 
second and either feet or yards per second, for 
Projectile Types 1, 2 and 8. (See References 7, 

8 and 9.) 

e . Form Factor 

It is not practical to determine the drag co- 
efficient as a function of Mach number for all 
projectiles. Instead, the drag coefficient of a pro- 
jectile whose shape differs only slightly from one 
of the typical projectiles is assumed to be pro- 
portional to the typical drag coefficient. If K Dt 
is the drag coefficient of a Typical Projectile, the 
form factor of a projectile whose drag coefficient 
is Kd is the ratio 

i t = K D /K Dt . (21) 

The form factor may be determined from the 
time of flight at short ranges, or from the range 
to impact or the position of burst at a few eleva- 
tions. If it is desired to estimate the form factor of 
a projectile without firing it, the Ordnance Corps 
Pamphlet, “ Design for Control of Flight Char- 
acteristics” 2 should be consulted. It explains the 
effects on drag of variations in length and radius 
of head, meplat diameter, base area, length of 
shell, and yaw. 

Generally, the form factor of a projectile with 
a square base and an ogive less than 1.75 calibers 
high should be referred to Projectile Type 1 ; that 
of a projectile with a square base and an ogive 
more than 1.75 calibers high, to Projectile Type 
8; that of a projectile with a boattail, to Projectile 
Type 2; and that of a finned projectile, to the 
T108 Shell. However, there are exceptions to 
these rules :e.g., K Dl is used for fin-stabilized mor- 
tar shells at low velocities, and Kdj fits the calcu- 
lated drag coefficient for some longer fin-stabilized 
shells. The drag coefficients of bombs are treated 
elsewhere. 10 

0.076,475 lb/ft* and the velocity of sound is 1116.89 fps; 
also the viscosity is 1.205 x 10“ * Ib/ft-sec, but the effect 
of Reynolds number is not considered in computing tra- 
jectories. 
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f. Ballistic Coefficient 

The ballistic coefficient relative to the tabulated 
drag coefficient is defined by the formula 

C t = m/itd’, (22) 

where m is the mass of the projectile in pounds, 
d is the caliber (or maximum diameter) 
in inches. 

g. Air Density 

The air density is a function of altitude. The 
exponential function adopted by the Ordnance 
Corps* is tabulated in References 7 and 11. The 
ratio of the density at an altitude y to that at the 
surface is H(y). 

4. Description of Typical Projectiles 

This pamphlet contains data that pertain spe- 
cifically to four types of projectiles: Types 1, 2 
and 8 and the 90-mm HEAT Shell T108. They 
can be applied to other shapes by using a form 
factor, as explained in paragraph 3e. 

a. Projectile Type 1 

Most of the service projectiles used by the 
European artillery during die last third of the 
nineteenth century had square bases and ogival 
heads of various apical angles. In Russia, during 
the years 1868 and 1869, General Magenski con- 
ducted some experiments at St. Petersburg to de- 
termine the law of air resistance of projectiles with 
a semi-apical angle of 48°22' (1.49-caliber radius 
head) at velocities from 172 to 409 m/s. In Eng- 
land, from 1866 to 1880, Bashforth conducted sim- 
ilar experiments with projectiles of the same form 
(it may be that the tables were reduced to this 
form) at velocities from 148 to 825 m/s. In Hol- 
land, Colonel Hojel conducted a few such experi- 
ments during 1883; the exact form is uncertain, - 
but the ogive was evidently somewhat longer than 
those of the Russian experiments. In Germany, 
the Krupp factory conducted a huge number of 
air resistance experiments at Meppen Proving 
Ground, starting in 1881 (die principal results 
were reported by Siacd in. the “Rivista” for 
March 1896) ; the velocities varied from 368 to 
910 m/s, and three different types of projectiles 

* The tabulated trajectory data based on the old stand- 
ard atmosphere give results tint are nearly die same as 
those based on the new atmosphere providing a suitable 
change is made in the ballistic coefficient. For example, a 
280-mm shell, fired with a muzzle velocity of 1795 fps, at 
an elevation of 45°, attains a range of 15790 yds.; the 
ballistic coefficient is 3.62 with ICAO atmosphere and 
3.50 with the old Ordnance Corps atmosphere. This func- 
tion was used in co m p u t in g trajectories before 1956, when 
the Department of Defense adopted the standard atmos- 
phere of the International Civil Aviation Organization. 


were used : one was the same as the Russian (or 
the results reduced to the same), but the others 
probably had longer ogives. In France, a com- 
mission of the naval artillery conducted many such 
experiments at Gavre Proving Ground from 1873 
to 1898 at velocities from 118 to 1163 m/s; they 
fired ogival projectiles with three semi-apical 
angles: 45°55' (1.64 CRH), 41°4 O' (1.98 CRH), 
and 31°45' (3.34 CRH). For further details, see 
References 12 to 16. 

In analyzing the results of the experiments, 17 
the Gavre Commission of Experiments assumed 
that the drag could be expressed in the form 

D = A a «v*f(V), (23) 

where 

D — is the drag (kg), 

A — the density of the air (kg/m*), 

g — the gravitational acceleration (m/sec*), 

a — the diameter of the largest cross section of 
the projectile (m), 

V — the velocity (m/s). 

The drag was' treated as a function of velocity 
without taking account of the velocity of sound, 
although the dependence of drag on the velocity of 
sound was strongly suspected: the coefficient f 
was plotted as a function of V alone. After 
studying the results obtained with projectiles of 
various ogival angles, the Commission confirmed 
the law deduced by Hdlie in 1888 that f(V) is 
proportional to the sine of die ogival angle, at 
least for values from 40° to 90°, within die experi- 
mental errors. 

In 1917, Chief Engineer Gamier of die Gtvre 
Commission published a table 1 * of log B and the 
corresponding logarithmic d e ri vativ e, with ve- 
locity as argument: this B is equivalent to f(V) 
for an ogival projectile with a 2-caliber radius 
head. This idol shape is what we call Projectile 
Type 1 (see Figure 1). The Ordnance Depart- 
ment of the United States Army then prepared a 
table of the logarithm of the drag f un ctio n 1 * 

G = kBV, ( 24 ) 

with V expressed in meters per second, and die 
factor k to take account of differences in defini- 
tions and standards. Below 600 m/s, k = 0.001 14 ; 
but at higher velocities, it is variable. G was 
called the Gavre function; later, die symbol was 
changed to Gi to denote that it pertains to Pro- 
jectile Type 1. 

In the table of log G, the argument, V*/100, 
varies from to 0 to 33,000 (equivalent to 1817 
m/s or 5961 fps). G has been tabulated for 
V*/100 up to 25,000 (1581 m/s or 5187 fps). 
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The product Gv has also been tabulated as a 
function of v in feet per second up to 6000 fps. 20 

b . Projectile Type 2 

A sketch of Projectile Type 2 is shown in 
Figure 1. This projectile has a 6° boattail y 2 
caliber long and an ogivo-conical head 2.7 cali- 
bers long. The experimental 4.7-inch Shell El 
had this shape. However, the resistance firings 
were actually conducted with two types of 3.3-inch 
shell: Type 155 with a 5° boattail and Type 157 
with a 7° boattail. The results were so close that 
they were grouped together to represent one type. 

These projectiles were fired from 1922 to 1925 
at velocities from 655 to 3187 fps. The results of 
these firings were fitted by the tabulated drag 
function G 2 (this is defined in paragraph 3d). 
Later, it was extrapolated to 10,000 fps by means 
of theoretical considerations. For the sake of 
greater smoothness and accuracy, two revisions 
have been made: G 2 . 1 , tabulated to 6000 fps, and 
G 2 . 2 , to 7000 fps. G 2.1 is approximately the same 
as G 2 , but G 2 .2 is lower at velocities above 2000 
fps. 

c. Projectile Type 8 

A sketch of Projectile Type 8 is shown in Fig- 
ure 1. This projectile has a square base and a 
secant ogive 2.18 calibers high. Its ogival radius 
is 10 calibers, which is twice that of a tangent 
ogive of the same height. 

The drag function G 8 was calculated from data 
obtained by the British Ordnance Board from 
three sources: 

( 1 ) Measurements of a model in the National 
Physical Laboratory wind tunnel at Mach numbers 
from 0.22 to 0.89. 

(2) Resistance firings of 2-pounder shell at 
Mach numbers from 0.65 to 2.8. 

(3) Resistance firings of 3-inch shell at Mach 
numbers from 0.35 to 4.15. 

Gg was first tabulated to 5000 fps. Later, a 
revised table, called G 8 .i, was extended to 6000 
fps. The corresponding drag coefficients are ap- 
proximately proportional to the British “1940 
Law," which was first tabulated to 5000 fps and 
later revised above 3750 fps and extended to 6000 
fps by a power law. 

d. HEAT Shell, 90-tnm, T108 

A sketch of the 90-mm High Explosive Anti- 
tank Shell T108 is shown in Figure 1. This pro- 
jectile has fins one caliber in span at the end of a 
long boom and a conical ogive about 2 calibers 
high. 

The drag of this projectile was determined at 
velocities from 700 to 2700 fps by free flight firings 


in the large Spark Range of the Exterior Bal- 
listic Laboratory at Aberdeen Proving Ground. 
The drag coefficient (defined in paragraph 3c) 
was fitted with a series of analytical expressions, 
which were used in tabulating it as a function of 
Mach number from 0.1 to 2.7. 

5. Trajectories in Air 

a , Equations of Motion 

Under standard conditions, the equations of mo- 
tion are 

GH ... GH . 

X = j~— x, y = g— y — g, (25) 

where x is the horizontal range, 
y — the altitude, 

H — the air density ratio, 
g — the gravitational acceleration (32.152 
ft /sec 2 ), 

and dots denote derivatives with respect to time. 
The initial conditions are 

Xo = 0, y c = 0, x©= v 0 cos * G , y<> = v c sin *& 0 , 

where v Q is the initial velocity, 

$ 0 — the angle of departure. 

These equations cannot be solved analytically, 
but they can be solved by a method of numerical 
integration explained in References 21, 22, 23 and 
24. 

b . Characteristics of Trajectories 

(1) Ballistic Tables, Trajectory data have 
been tabulated for Projectile Types 1 and 2 and 
for some bomb shapes. 

The “Exterior Ballistic Tables Based on Nu- 
merical Integration" (References 7, 25, and 26) 
pertain to Projectile Type 1. The meter is the 
unit of length. The trajectories were computed 
forward and backward from the summit and cut 
off at several altitudes. If C. is the “summital 
ballistic coefficient" of a trajectory, and y. the max- 
imum ordinate of a partial trajectory, the ballistic 
coefficient C of this partial trajectory satisfies the 
relation 

logxoC = log 10 C g — 0.000,045 y g . (26) 

For each y., and summital velocity v„ 
Volume I gives the horizontal range, time, and 
velocity components at the beginning and end 
of the partial trajectory. If the muzzle velocity, 
angle of departure, and ballistic coefficient are 
known. Volume II may be used to find the 
time of flight to the summit, t*, the velocity at the 
summit, the horizontal range to the summit, x„, 
and the maximum ordinate. Then C, can be com- 
puted and trajectory data can be found in Volume 
I for the given C. and v. at various altitudes. 
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puted and trajectory data can be found in Volume 
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However, terminal values can be found more con- 
veniently from Volume III, which gives the hori- 
zontal range X, time of flight T, angle of fall 
<*>, and terminal velocity v w ; the arguments are 
the same as in Volume II : the muzzle velocity 
varies from 80 to 920 meters per second, the 
angle of departure from 0° to 80° (at intervals of 
1° from 0° to 30°, and 2° from 30° to 80°), the 
ballistic coefficient from 1.000 to 17.78 (at inter- 
vals of 0.05 in log C). 

The “Brief Exterior Ballistic Tables for Pro- 
jectile Type 1” (Reference 27) contains values 
that were obtained by interpolation in Volume 
III and conversion to English units (yards and 
feet per second ) . The muzzle velocity varies from 
200 to 3000 feet per second in 100 foot per second 
increments. The angle of departure varies from 
15° to 75° in 15° increments. Obviously, inter- 
polation for angle of departure in this table is 
difficult ; nevertheless, the table is useful for many 


purposes. The table of range for an elevation of 
45°, which is nearly the maximum range, is repro- 
duced here (Table I). 

The “Exterior Ballistic Table for Projectile 
Type 2 99 (Reference 28) consists of three parts: 
Table I, Summital and Terminal Values; Table 
II, Coordinates as Functions of Time of Flight; 
and a large graph of Maximum Range vs. Ballistic 
Coefficient. The tables have been extended to 
include data for every 15° from 0° to 90°. The 
muzzle velocity varies from 1200 to 5200 feet per 
second, and the ballistic coefficient from 0.50 to 
10.00. The coordinates are given in yards, the 
time in seconds, and the velocity in feet per second. 
Table I gives the time to the summit, the hori- 
zontal range to the summit, the maximum 
ordinate, the time of flight, the total range, the 
angle of fall, and the terminal velocity. Table II 
gives the horizontal range and altitude at every 
5 seconds of time. Unfortunately, the velocity 


TABLE I 

MAXIMUM RANGE FOR PROJECTILE TYPE 1 
ELEVATION 45° 
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404 

405 

406 

407 

408 

408 

409 

409 

300 

847 

856 

864 

871 

877 

883 

888 

892 

897 

900 

904 

907 

<uo 

■EZE1 

400 

1419 

1441 

1462 

1481 

1499 

1514 

1528 

1541 

1553 

1564 

1573 

1582 

1590 

•I KM 

500 

2077 

2121 

2162 

2201 

2237 

2270 

2300 

2328 

2354 

2378 

2399 

2418 

2436 


(500 

2786 

2862 

2934 

3003 

3067 

3126 

3181 

3232 

3279 

3322 

3361 

3397 

3428 

3461 

700 

3521 

3(538 

3749 

3856 

3958 

4053 

4140 

4222 

4301 

4372 

4438 

4500 

4556 

4606 

800 

4247 

4411 

4570 

4724 

4870 

5008 

5141 

5267 

5384 

5494 

5595 

5690 

5777 

5859 

900 

4920 

5136 

5349 

•5556 

5755 

5947 

6133 

6310 

6476 

6634 

6783 

6923 

7054 

7174 

1000 

5489 

5757 

6024 

6286 

6542 

6792 

7037 

7272 

7498 

7714 

7919 

8115 

8299 

8472 


1RE1M 

6238 

6553 

6866 

7176 

7482 

7784 

8079 

8366 

8644 

8913 

9172 

9420 

9658 




6958 

7314 

7671 

8027 

8381 

8733 

9077 

9417 

9749 

10073 

10388 

10691 




7293 

7689 

8087 

8488 

8891 

9292 

9693 


10489 

10872 

11253 

11626 


6759 

7172 

7594 

8025 

8460 

8905 

9353 

9804 

10258 


11165 

11616 

12065 

12508 


6979 

7419 

7872 

8336 

8810 

9294 

9787 

10287 

10792 

11301 

11814 

12329 

12845 

13360 

U500 

7185 

7651 

8134 

8630 

9141 

9664 


10747 

11302 

11869 

12442 

13021 

13605 

14194 

1700 

7381 

7873 

8384 

8912 

9458 


10599 

11193 

11798 

12421 

13056 

13701 

14356 

15019 

1800 

7569 

8088 

8627 

9187 

9767 


10989 

11630 

12288 

12d67 

13662 

14374 

15099 

15839 

1900 

7752 

8296 

8863 

9454 



11371 

12059 

12771 

13505 

14262 

15041 

15839 

16659 

2000 

7931 

8498 

9093 

9715 



11749 

12483 

13247 

14037 

14859 

15706 

16580 

17480 

2100 

8106 

8698 


9972 

10655 

11372 

12121 

12903 

13719 

14567 

15454 

16372 

17323 

18306 

2200 

8278 

8894 

9542 

10225 

10942 

11697 

12489 

13319 

14188 

15096 

16049 

17038 

18068 

19136 

2300 

8447 

9086 


10473 

11225 

■ME1 

12853 

13732 

14655 

15624 

16644 

17706 

18820 

19975 

21 ■' 0 

8612 

9274 

9975 

10718 

11.505 

12336 

13215 

14143 

15121 

16152 

17239 

18375 

19574 

20821 

2500 

8774 

9459 

10188 

10961 

11782 

12652 

13675 

14552 

15585 

16679 

17834 

19048 

20330 

21675 

2600 

8935 

9642 

10397 

11201 


12965 

13931 

14958 

16048 

17205 

18430 

19725 

21095 

22539 

2700 

9092 

9823 

10603 

11438 

13237 

13275 

14285 

15362 

16510 

17731 

19027 

20405 

21866 

23412 

2800 

9247 

10001 

10807 

11672 

12595 

13582 

14637 

15764 

16970 

18257 

19625 

21091 

22645 

24294 

2900 

9399 

10177 

11010 

11903 

12861 

13885 

14987 

16165 

17429 

18782 

20225 

21778 

23428 

25184 

3000 

9549 

10351 

11210 

12132 

13123 

14187 

15334 

16563 

17885 

19303 

20823 

22464 

24213 

mu 
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TABLE I (Continued) 

MAXIMUM RANGE FOR PROJECTILE TYPE 1 
ELEVATION 45° 


Range in yards. Muzzle Velocity (v 0 ) in fps 


{§ 




0.750 

5.623 

0.800 

6.310 



0.950 

8.913 

1.000 

10.00 



1.150 

14.13 

1.200 

15.85 

1.250 

17.78 

HI 

409 


409 

m 

410 

410 

411 

411 

412 

412 

412 

413 

413 

413 



913 

915 

■ I 

918 

919 

921 

923 

924 

925 

925 

926 

927 

928 

1“ 


1597 

1604 


1614 

1619 

1624 

1627 

1630 

1634 

1637 

1639 

1641 

1643 

Eaufl 

2436 

2452 

2467 

2480 

2491 

2501 

2511 

2519 

2526 

2533 

2540 

2545 

2550 

2555 

600 

3428 

3461 

3489 

3513 

3535 

3555 

3575 

3591 

3605 

3619 

3631 

3640 

3649 

3659 

700 

4556 

4606 

4651 

4695 

4733 

4766 

4800 

4829 

4855 

4879 

4899 

4918 

4934 

4951 

800 

5777 

5859 

5934 

6001 

6064 

6120 

6172 

6219 

6260 

6298 

6333 

6364 

6394 

6418 

900 

HmI 

7174 

7284 

7389 

7484 

7573 

7652 

7725 

7791 

7852 

mriOfM 

7958 

miwm 

8044 

1000 

8299 

8472 

8634 

8786 

8928 

9060 

9180 

9292 

9394 

9489 

9575 


9724 

9787 

1100 


9658 

9882 

10094 

10295 

10484 

10660 

10825 

10979 


11252 

11373 

11485 

11588 

1200 



10984 

11266 

11536 

11794 

12038 

12268 

12487 

12691 

12882. 


13228 

13382 

1300 . 

11253 

11626 

11990 

12344 

12688 

13020 

13339 

13644 

13936 

14213 

14476 

14723 

14958 

15175 

1400 


12508 

12945 

13375 

13796 

14207 

14604 

14991 

15364 

15722 

EEX5I 

n6394 


■E*a 

1500 

12845 


13872 

14381 

14881 

15374 

15857 

16331 

16790 

17235 

17667 


18474 

18852 

1600 


14194 

14783 

15370 

15956 

16536 

17109 

17674 

18227 

18766 

19291 

19798 


20756 

1700 

14356 


15686 

KOMI 


17702 

18371 


19682 


20949 

21559 



1800 

wmm 

15839 

16588 

17344 


18875 

19643 



21913 





1900 

15839 

16659 

17491 

18338 

19198 

20063 

20936 


22678 

wzm l 






16580 


18399 



21270 

22253 

23243 

24222 

25222 





2100 

17323 

18306 

19315 


21415 

22497 


24710 

25829 

26057 





2200 


19136 

20240 


22551 

23748 

24972 

26217 







2300 

18820 

19975 

21177 

22421 

mssa 

msMm 


27771 







2400 

19574 

20821 

22124 

23478 

24883 

26334 

27832 

29372 







2500 

20330 

21675 

23082 

24553 


27673 

29321 

31028 







\mm 

21095 

22539 

24056 

25646 

27307 



32748 







2700 

21866 

23412 

25045 

26763 

28564 


32446 









22645 

24294 

26040 


29849 

KtTTTjj 

1 mimm 









23428 

25184 

27049 

29037 

31183 


35797 









24213 

26083 

28080 



34944 

37558 









TABLE II 

MAXIMUM RANGE OF’PROJECTILE TYPE 2 
Range in Units of 10 Yards 


Muzase Velocity 
fpa 

Log C, 

-0.3 

-0.2 

-0.1 

0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

1200 




709 

861 

1004 

1129 

1227 

1308 

1600 


.... 

♦ • ♦ . 

851 

1088 

1335 

1584 

1821 

2030 

2000 




1001 

1322 

1690 

2091 

2409 


2400 




1143 

1561 

2071 

2659 

3295 

3983 

2800 

759 

905 

1078 

1284 

1806 

2479 

3299 

4322 

5359 

3200 

819 

983 

1182 

14% 

2062 

2926 

4074 

5589 

7006 

3600 

877 

1061 

1287 

1571 

2335 

3438 

5061 

7201 

9034 

4000 

934 

1139 

1394 

1720 

2629 

4056 

6329 

9106 

11356 

4400 

992 

1218 

1504 

1874 

2959 

4834 

7917 

11335 

14019 

4800 

1051 

1299 

1617 

2037 

3342 

5845 

9836 

13888 

17015 

5200 

i 

1112 

1384 

1735 

2220 

3795 

7196 

12091 

16752 

20333 


components are not tabulated. The graph has a 
curve of maximum range vs. ballistic coefficient 
for every 400 feet per second muzzle velocity; it 
also lists the angle of departure for maximum 
range. A table for the maximum range is included 


in this pamphlet (Table II). 

Some trajectories have been computed for Pro- 
jectile Type 8 at an angle of departure of 45°. 
The range, which is nearly the maximum range, 
is tabulated here (Table III). 
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TABLE III 

MAXIMUM RANGE FOR PROJECTILE TYPE 8 
ELEVATION 45° 

Range in Units of 10 Yards 


Muzzle Velocity ips 

Log C 8 

-0.2 

q 

© 

0.2 

0.4 

0.6 

0.8 

1.0 

1200 

521 

665 

815 

961 


1201 

1288 

WIT? 1 

680 

923 

1223 

1581 

1934 

2403 

2826 

mmwmmmm 

822 

1168 

1650 

2287 

3124 

4231 

5251 


964 

1432 

2145 

3162 

4865 

6936 

8741 

4400 

1108 

1720 

2726 

4470 

7492 

10884 

13661 

5200 

1258 

2040 

3424 

6604 

11693 

16582 

20319 


When bombs were dropped at low speeds, tra- 
jectory data based on the Gavre drag function Gi 
were accurate enough for calibrating the sights. 
In fact, the descending branches of the trajectories 
in Volume I of the “Exterior Ballistic Tables ,, 
could be used for this purpose. Actually, some 
“Bomb Ballistic Reduction Tables” were pub- 
lished, which were more convenient for this pur- 
pose. However, at transonic speeds, the drag 
functions of modern bombs differ considerably 
from Gi. Drag coefficients have been determined 
for several shapes, but three of them have been 
selected as bases of “Bomb Ballistic Tables” (Ref- 
erence 10) : drag coefficient A was determined 
from wind tunnel tests of a shape designed by the 
National Advisory Committee for Aeronautics, and 
is probably the lowest drag coefficient that any 
bomb will have ; drag coefficient B was determined 
from range bombing of the 3000-pound Demoli- 
tion Bomb Ml 19, and is suitable for other mem- 
bers of the series including the 750-pound Bomb 
Ml 18 and the 10,000-pound Bomb M121; drag 
coefficient C was determined from resistance fir- 
ings of the 75-mm Slug Mark I, and is probably 
the highest coefficient that any bomb will have. 

The arguments of the “Bomb Ballistic Tables” 
are release altitude ( 1000 to 80,000 feet above sea 
level), true air speed of release in level flight (250 
to 2000 knots), and performance parameter (0 
to 1.5). The performance parameter P is inverse- 
ly proportional to the ballistic coefficient ; it is de- 
fined by the formula 

P = 100d 2 /m, (27) 

where d is the maximum diameter of the bomb in 
feet. 


m — the mass of the bomb in pounds. 

The tables give the horizontal range from release 
to impact (feet), the time from release to impact 
(seconds), the striking velocity (ft/sec), and the 
striking angle (degrees). The drag coefficients 
are tabulated as functions of the Mach number in 
the appendix ; they are graphed in the intro- 
duction. 

(2) Graphs . Figures 3, 4 and 5 are graphs of 
the total time of flight for an elevation of 45°. The 
time of flight is useful in estimating wind effects. 

For antiaircraft firing and for firing to the 
ground at different altitudes, points along the tra- 
jectory are needed. Figures 6 and 7 show the 
trajectories computed with G 2 and Gg at angles of 
departure of 45° and 60° (there is a different 
graph for each muzzle velocity, each elevation, 
and each Projectile Type). Trajectories based 
on Gi are not included because no antiaircraft pro- 
jectiles have a blunt head. The T108 HEAT 
Shell is usually used only in direct fire at short 
ranges. Trajectory data for other elevations can 
be estimated from those at 45° or 60° by assuming 
that the slant range is the same for a given time of 
flight. 

Unlike the vacuum trajectories, in air, the 
horizontal range to the summit is more than half 
the total range, and the time to the summit is lesjs 
than half the total time; however, the maximup 
ordinate is approximately (g/8)T 2 . In air, the 
terminal velocity is less than the initial velocity 
and the angle of > fall is greater than the angle of 
departure. At a given elevation, the total r^nge 
increases less rapidly than the square of the 
muzzle velocity; at a given muzzle velocity, the 
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range is a maximum at an angle somewhat differ- 
ent from 45°: lower for small ballistic coefficients, 
higher for large ones. However, the range at 45° 
is nearly as much as the maximum. 

6. Siacci Tables 
a. Spin-stabilized Projectiles 

(1) Tabular Functions . In 1880, the Italian 
Col. Siacci introduced the Space, Time, Altitude, 
and Inclination Functions, (Reference 29), which 
simplify the calculation of trajectory data under 
certain restrictions. In particular, their use will 
be explained for three conditions: ground fire at 
low superelevations, antiaircraft fire (where the 
superelevation is small), and aircraft fire. These 
and some auxiliary functions based on Gi, G 2.2 
(the second revision of G 2 ) and Gg.i (the first 
revision of Gg) have been tabulated for the Air 
Force (References 30, 31 and 32). Since these 
Siacci tables were arranged primarily for use in 
computing aircraft firing tables, they were called 
Aircraft Tables; however, they can be used for 
other purposes, as explained below, in lieu of other 
forms. In these Aircraft Tables the argument is 
the component u along the line of departure of the 
velocity relative to the air, which is assumed npo- 
tionless. The following functions and some of 
their derivatives with respect to u are tabulated: 


Space 

o- 1 

r u 

du 

(28) 

p -J 


G(u) ’ 

Time 

*=J 

r u 

du 

(29) 


uG(u) ’ 

Inclination* 

q ‘=J 

r u 

g du 

(30) 

L 

u 2 G(u) * 


Altitude* q = f rt \ 9 

J u G(u) 

(31) 

Velocity increment A'u = — ^G(u), 

(32) 

(K D 5 = 16.4 per rad 2 ), 


Damping coefficient c" = G(u) /2u. 

(33) 

The upper limit of integration, U, is 6000 fps for 
Projectile Type 1 and 7000 fps for Types 2 and 


8. K is the yaw-drag coefficient, defined by 
the formula 

Kd = Kd^ ( 1 + K d S 2 ), (34) 

where is the drag coefficient for 0 yaw, 

8 — is the yaw in radians. 

(2) Ground Fire at Low Elevations 

(a) Formulas . Although the Aircraft Tables 
are not accurate in computing complete trajectories 
at high elevations, they can be used to obtain tra- 
jectory data at superelevations below 15° (for the 
derivation and application of the following for- 
mulas, see Reference 33). Here, we shall take 

p/p o — as the ratio of air density to normal, as- 
sumed constant for each problem, 

a/a© — the ration of the velocity of sound to 
normal, assumed constant for each 
problem. 

If F is the temperature in degrees Fahrenheit, 
a/ao = V(459 + F)/518. (35) 

Then the horizontal range is 

(36) 


X = (Cpo/p) cos 0. [p(ua«/a) — p(uoa«/a)]. 


the time of flight is 

t,= (Cpoao/pa) [t(uao/a) — t(tM./a)], 

the altitude is 

y = xtan — (Cp«a./pa) a [q(ua«/a) — q(n^/&)] + x(Cp.ao/pa) sec*. [qi(u 0 a./a)], 
and the slope of the trajeetoiy is 

tan * = tan *. — (Cpoa.*/pa 4 ) sec *» [qi (ua«/a) — qi (uoa<,/a)]. 


The angle of departure for y = 0, may be found by formula 


j _ Cppap r q(ua./a) — q(uoao/a) 
0 pa Lp(ua«/a) — p(u«a./a) 


qi(u«ao/a)J. 


If p/p„ = 1 and a/a. = 1, let p = p(u), p. = p (u.), etc. 


(37) 

(38) 

(39) 

(40) 


* In most places, the altitude and inclination functions are defined as twice these integrals. 
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range is a maximum at an angle somewhat differ- 
ent from 45°: lower for small ballistic coefficients, 
higher for large ones. However, the range at 45° 
is nearly as much as the maximum. 

6. Siacci Tables 
a. Spin-stabilized Projectiles 

(1) Tabular Functions . In 1880, the Italian 
Col. Siacci introduced the Space, Time, Altitude, 
and Inclination Functions, (Reference 29), which 
simplify the calculation of trajectory data under 
certain restrictions. In particular, their use will 
be explained for three conditions: ground fire at 
low superelevations, antiaircraft fire (where the 
superelevation is small), and aircraft fire. These 
and some auxiliary functions based on Gi, G 2.2 
(the second revision of G 2 ) and Gg.i (the first 
revision of Gg) have been tabulated for the Air 
Force (References 30, 31 and 32). Since these 
Siacci tables were arranged primarily for use in 
computing aircraft firing tables, they were called 
Aircraft Tables; however, they can be used for 
other purposes, as explained below, in lieu of other 
forms. In these Aircraft Tables the argument is 
the component u along the line of departure of the 
velocity relative to the air, which is assumed npo- 
tionless. The following functions and some of 
their derivatives with respect to u are tabulated: 


Space 

o- 1 

r u 

du 

(28) 

p -J 


G(u) ’ 

Time 

*=J 

r u 

du 

(29) 


uG(u) ’ 

Inclination* 

q ‘=J 

r u 

g du 

(30) 

L 

u 2 G(u) * 


Altitude* q = f rt \ 9 

J u G(u) 

(31) 

Velocity increment A'u = — ^G(u), 

(32) 

(K D 5 = 16.4 per rad 2 ), 


Damping coefficient c" = G(u) /2u. 

(33) 

The upper limit of integration, U, is 6000 fps for 
Projectile Type 1 and 7000 fps for Types 2 and 


8. K is the yaw-drag coefficient, defined by 
the formula 

Kd = Kd^ ( 1 + K d S 2 ), (34) 

where is the drag coefficient for 0 yaw, 

8 — is the yaw in radians. 

(2) Ground Fire at Low Elevations 

(a) Formulas . Although the Aircraft Tables 
are not accurate in computing complete trajectories 
at high elevations, they can be used to obtain tra- 
jectory data at superelevations below 15° (for the 
derivation and application of the following for- 
mulas, see Reference 33). Here, we shall take 

p/p o — as the ratio of air density to normal, as- 
sumed constant for each problem, 

a/a© — the ration of the velocity of sound to 
normal, assumed constant for each 
problem. 

If F is the temperature in degrees Fahrenheit, 
a/ao = V(459 + F)/518. (35) 

Then the horizontal range is 

(36) 


X = (Cpo/p) cos 0. [p(ua«/a) — p(uoa«/a)]. 


the time of flight is 

t,= (Cpoao/pa) [t(uao/a) — t(tM./a)], 

the altitude is 

y = xtan — (Cp«a./pa) a [q(ua«/a) — q(n^/&)] + x(Cp.ao/pa) sec*. [qi(u 0 a./a)], 
and the slope of the trajeetoiy is 

tan * = tan *. — (Cpoa.*/pa 4 ) sec *» [qi (ua«/a) — qi (uoa<,/a)]. 


The angle of departure for y = 0, may be found by formula 


j _ Cppap r q(ua./a) — q(uoao/a) 
0 pa Lp(ua«/a) — p(u«a./a) 


qi(u«ao/a)J. 


If p/p„ = 1 and a/a. = 1, let p = p(u), p. = p (u.), etc. 


(37) 

(38) 

(39) 

(40) 


* In most places, the altitude and inclination functions are defined as twice these integrals. 
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Then these formulas can be expressed more 


simply : 

x = C cos (p — Po) , (36a) 

t f = C(t — 1 0 ), (37a) 

y ~ x tan $ 0 — C 2 (q — q 0 ) + 

C sec *S 0 x qio) (38a) 

tan $ = tan # 0 — C sec # 0 (qi — qio) f (39a) 
and, for y = 0, 

sin«. = C (40a) 

LP Po J 


(b) Example 

Data : Tank Gun, 
M339. 

Caliber 
Mass 

Muzzle Velocity 
Form Factor 
Air Density Ratio 
Sound Velocity Ratio 


M1A2; Shot, AP, 

3.000 in. 

14.58 lb. 

3200 fps 
1.164 on G«.t 

1.000 
1.000 


76-mm, 


Problem : To find the range, time of flight, angle 
of departure, and angle of fall where 
the velocity is 2000 fps. 

Solution : From the Aircraft Table Based on G $.\ : 


u 

P 

t 

q 

qi 

2000 

23079.3 

6.5132 

489.0 

.06523 

3200 

15425.1 

3.4997 

161.6 

.02637 

Diff. 

7654.2 

3.0135 

327.4 

.03886 


C 8 .i = 14.58/1.164 (3.000) 2 = 1.392 (Eq. 22) 

sin 4 0 = 1.392 (327.4/7654.2 - .02637) = 
.02283 (Eq. 40a) 

•So = 23.26 mils, cos *S 0 = .99974, sec. -S 0 = 
1.00027 

tan -S 0 = .02284 

1.392(1.00027) (.03886) = .05411 

tan *& = —.03127 (Eq. 39a) 

•S — — 31.84 mils 

x = 1.392 (.99974) (7654.2) = 

10652 ft (Eq. 36a) 

t f = 1.392 (3.0135) = 4.19 sec (Eq. 37a) 


Note: Since u is the component of the velocity 
along the line of fire, the tangential velocity is 
v = u cos $ 0 /cos £. 

In this example, v = 2000 (. 99974 )/.99951 = 
2000.5 fps. 

(3) Antiaircraft Fire 

(a) Formulas . These tables can also be 
used for the flat part of antiaircraft trajectories, 
where the superelevation is less than 15° (for the 
derivation and application of the following formu- 
las, see Reference 33). Let p/p 0 and possibly C 
be known variable functions of x, but assume a 
constant value of a/a* (either 1 or an average). 
Since p/p 0 is a function of y, denoted by H, its 
relation to x can be found by the approximate sub- 
stitution 

y *=*=* x tan # 0 . (41) 

Let m = tan m 0 = tan •&<>. 

Then the following trajectory data can be found 
by Simpson's rule or some other method of nu- 
merical integration: 


x 


p (uao/a) = p(uoao/a) + sec # 0 



dx, (42) 


o 


X 

t f = sec £ 0 j* “ , (43) 

o 


m 


g sec 2 


f g sec 

— m ° -J— 


dx, 


(44) 


X 


-f 


m dx. 


o 


( 45 ) 


Check : x tan = 10652 (.02284) = 243.29 

- Cla (q - q.) = (1-392)2 (327A) = - 634.40 

Qu sec \ X qio = 1.392 (1.00027) (10652) (.02637) = 391.11 

y = 


10 


0.00 (Eq. 38a) 



(b) Example 


1 Data: Automatic (AA) Gun, 40-mm, M2A1; Shot, AP, M81 







Muzzle Velocity 

2870 fps 









Ballistic Coefficient 

0.574 on Gm 








Angle of Departure 

1200 mils 









Air Density Ratio 

K 

1 

ffl 









Sound Velocity Ratio: 1.000 









Problem: To find the time of flight and altitude for horizontal ranges from 0 to 2400 ft. 




Solution: (See remarks below and use the Aircraft Table Based on G*.i). 





X 

Mhy 

H 

f Hdx 

p(u) 

u 

2613.1 


/2613.1V 

r/26i3.iy dx 

m 

fm 

ft 

eat 





u 

sec 

\ u ) 

J V U / 


ft 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


.00000 

.00903 

.01987 

1.0000 

.9774 

.9553 

0.00 

17439.6 

2870.0 

0.9105 

0.000 

0.8290 

0 

2.4142 

0 

586.49 

20109.5 

2446.4 

1.0681 

0.594 

1.1408 

59i 

2.3952 

1443 

n«»y« 

.02980 

.03973 

.9337 

.9126 

1146.76 

22660.1 

2060.6 

1.2681 

1.295 

1.6081 

im 

2.3687 

2872 

1500 

1800 

.04966 

.05960 

.8919 

.8718 

1681.96 

25096.6 

m6.3 

i.5225 

2.i32 

2.3180 

2593 

2.3308 

4282 

2100 

2400 

.06953 

.07946 

.8521 

.8328 

2193.26 

27424.2 

1413.7 

1.8484 

3.143 

3.4166 

43i4 

2.2755 

5664 


Remarks : 

Col 1. The interval of x should be chosen small 
enough so that the integrations may be performed 
accurately by Simpson's rule : not more than 600 
feet. 

Col 2. M = logio e s= 0.43429. h = 0.000,031,- 
58 ft -1 . The estimate of y = x tan & 0 . A 
slightly smaller factor would be better at the 
longer ranges. 

Col 3. For convenience, the formula H = anti- 
cologio Mhy was used. The density function 
H(y) has been tabulated (Reference 11) ; it gives 
the same values. The ICAO “Standard Atmos- 
phere — Tables and Data” 1 gives slightly different 
values. 

Col 4. Simpson's rule was used, but the trap- 
ezoid rule gives 2193.30 for the last value. 

Col. 5. The first value was found in the Aircraft 
Table for Uo — 2870 fps. The other values were 
computed by Eq. (42) with 

sec « 0 /Ci.i = 2.6131/0.574 = 4.5524. 

Col 6. The values of u after the first were found 
in the Aircraft Table by inverse interpolation, 
using the first two terms of Taylor’s series. 

Col 7. sec $ 0 = 2.6131. This is multiplied by 
1000 to avoid writing zeros. 


Col 8 . Eq. (43). Simpson's rule gives 3.131 
for the last Value, but the trapezoid rule was used 
in order to obtain intermediate values. 

Col 9 . The square of Col. 7. 

Col 10. Note the factor 10 6 in these values. 
Simpson's rule gives 4259 for the last value, but 
the trapezoid rule was used so as to obtain more 
values. 

Col 11. Computed by Eq. (44) with 

m 0 = tan £ 0 , g = 32.152 ft/sec 2 . 

Col 12. Computed by Eq. (45). Simpson's rule 
gives 5666 ft for the final value. 

(4) Aircraft Fire 

(a) Formulas. These tables are especially 
useful for aircraft fire (for the derivation and ap- 
plication of the following formuas, see References 
34 and 35). Before using them, the initial con- 
ditions must be determined. Let 

A — be the azimuth, measured from the direc- 
tion of motion of the airplane, 

Z — the zenith angle, measured from the upward 
vertical, 

v 0 — the muzzle velocity (relative to the gun), 
w — the true air speed, 

Uo — the initial velocity relative to the air, 

S 0 — the initial yaw relative to the motion of the 
air. 
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The initial velocity and yaw may be found by the 


formulas 

14 2 ss v 0 2 + w 2 -f 2wv 0 sin Z cos A, (46) 
So 2 = (1 — sin 2 Z cos 2 A) w 2 /tto 2 (rad) 2 . (47) 
In forward firing, A = 0 ; Z = 1600 mils, 

Uo = v 0 + w, (46a) 

8 0 = 0. (47a) 

In rearward firing, A = 3200 mils, Z = 1600 
mils, 

u 0 = v 0 — w, (46b) 

S 0 = 0. (47b) 

The air density ratio H and the sound velocity 


ratio a/ao may be taken as constants. Approxi- 
mately, if Y is the altitude of the airplane above 
sea level, 

H e“ hY , h = 0.000,031,58 ft" 1 , (48) 

a/ao = e“ kY , k = 0.000,003,01 ft" 1 . (49) 

If s* is the standard stability factor, which cor- 
responds to the muzzle velocity, muzzle spin, and 
standard surface air density, the initial stability 
factor is 

So = s,v 0 2 /Uo 2 H. ( 50) 


If h is the yawing moment damping factor, 
x — the cross wind force damping factor, 

c'=(h + x)/2u, (51) 


the loss in velocity due to yaw is 
At, — *° 3a (so — 1/2) A'u (up/a) 

a. (so— 1) Cc' + c" (u<#«/a) 
Then the effective initial velocity is 


(52) 


14 = Uo — Au. 


(53) 


Now we can find the distance measured along 
the line of departure to a point directly above the 
projectile, 

P = (C/H) [p(ua«/a) — p (u^/a )], (54) 

the drop of the projectile, 

Q = (G*o/Ha) 2 [q(uao/a) — q(u*ao/a)] 

— ((X/Ha) Pq x (u^/a), (55) 

and the time of flight is, 
t* = (Cao/Ha) [t(ua«/a) — t (u^/a)]. (56) 

The horizontal and vertical components, if de- 
sired, may be found by the formulas 
x = P cos -So, (57) 

y = P sin -S 0 — Q, (58) 

with the angle of departure determined by the 
formula 


tfln 2A — ~ 

° (sin Z sin A) 2 + (sin Z cos A -f- w/v 0 ) 2 
or 

cot 2 4 0 ss (tan Z sin A) 2 -f- (tan Z cos A + w/v 0 cos Z) 2 . 


(59) 

(59a) 


Illustrative example is given on the following page. 
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(b) Example 
Data: 

Azimuth : 1200 mils 

Zenith Angle : 800 mils 

Muzzle Velocity : 3000 fps 

Air Speed : 1000 fps 

Altitude : 30,000 ft 

Standard Stability Factor : 3.66 

Damping Factors h + * : 10.6 sec -1 

at 3000 fps 

Ballistic Coefficient : 0.235 on Gg,x 

Problem : To find the time of flight and drop for 
a Siacci range (P) at 1200 ft. 

Solution: (Use the Aircraft Table Based on 

Gg.i). 

cos A = 0.38268 


sin Z = 0.70711 
sin Z cos A ~ 0.27060 
2w v 0 = 2 x 1000 x 3000 = 6,000,000 
v 0 2 = (3000) 2 = 9,000,000 

w 2 = (1000) 2 = 1,000,000 

2wv 0 sin Z cos A = 1 ,623,600 

tio 2 = 1 1,623,600 (Eq. 46) 

Uo = 3409.3 fps 

5 0 2 = 0.72940 x 1,000,000/11,623,600 = 
0.06275 (Eq. 47) 

= 0.2505 rad (255.2 mils), (1 rad = 

1018.6 mils). 

hY =0.000,031,58 x 30,000 = 0.9474 
H = 0.3877 (Eq. 48) 
kY = 0.000,003,01 x 30,000 = 0.0903 
a/a. = 0.9137 (Eq. 49) 

So = 3.66 x 9,000,000/11,623,600 x 0.3877 = 
7.31 (Eq. 50) 

c' = 10.6/2 x 3000 = 0.00177 ft- 1 
u^/a = 3409.3/0.9137 = 3731.3 fps 

A'u(uoao/a) = 1.4003 sec -1 
(from Aircraft Table) 

c"(u<>ao/a) = 2288 x 10“® ft” 1 
(from Aircraft Table) 


u. = 3409.3 - 206.5 = 3202.8 fps (Eq. 53) 
Cga/H = 0.235/0.3877 = 0.6061 
CgWHa = 0.6061/0.9137 = 0.6633 
(Cg.iao/Ha) 2 = (0.6633)* = 0.4400 
PH/Cs.1 = 1200/0.6061 = 1979.9 ft 
u*a«/a = 3202.8/0.9137 = 3505.3 fps 
Find p(u«a«/a) in Table, then p(uao/a) by Eq. 
54. 


uao/a 

P 

t 

q 

q« 

3174.3 

15580.5 

3.548 

165.7 

.02686 

3505.3 

13600.6 

2.955 

118.4 

.02114 


1979.9 

0.593 

47.3 

.00572 


t f = 0.6633 x 0.593 = 0.393 sec (Eq. 56) 

0.4400 x 47.3 = 20.8 

0.6633 x 1200 x 0.02114 = 16.8 
Q = “T0ft (Eq. 55) 

(c) Interpolation . The derivatives listed it 
the Aircraft Tables should be used for inter- 
polating backward or forward from the nearest 
tabular value of the argument. For instance, in 
the last example, it is desired to find p for u = 
3505.3 and then to find u for p = 15580.5. First, 
the table gives 


at u = 3510, 



Hence, Au = — 4.7, 

and, at u = 3505.3, 

Then we have p = 15580.5 
and the table gives 

p = 15606.2, - 6.044 

du 

Hence, Ap = — 25.7 

and 


p = 13572.81, 


Ap — + 2/. 84, 
p = 13600.65. 


at u = 3170. 

Au = +4.3, 
u = 3174.3. 


Au ~ 0.06275 x 0.9137 ^ x 0 235 x 0 .00177 + 0.000,022,88 
= M5733 ^ (E * 52) 
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b„ Fin-stabilised Projectiles 

(1) Tabular Functions . Reference 6 contains 
a table of modified Siacci functions based on the 
drag coefficient for the 90-mm High Explosive 
Antitank Shell T108, which has rigid fins. These 
tables can be used for other similar projectiles, 
although they are not suitable for all fin-stabilized 
projectiles. The following functions and their 
derivatives are tabulated as functions of Mach 
number M: 


(2) Application 

(a) Formulas . The following formulas are 
applicable for computing trajectory data at super- 
elevations below 15°. 

Let 

m be the mass of the projectile in pounds, 
i — the form factor of the projectile, 

d — the maximum diameter of the projectile in 
feet, 


2.7 

Space S= f* MKd(M) . (60) 

M 

2.7 

Time T = /m€(M)’ (61) 

M 

2.7 

Inclination 1= j- f - . (62) 

M 

Al,i, “ d ' A = f-m % • <«> 

M 

The factor in the inclination function is 
g/ao = 32.1S4/1120.27 = 0.0287 sec” 1 . 


v — the velocity of the projectile in feet per 
second, 

v 0 — the muzzle velocity in feet per second, 
a — the velocity of sound in air in feet per second, 
p — the density of air in pounds per cubic foot, 
# — the angle of inclination of the trajectory, 

*o — the angle of departure. 

Also let 

a G — 1120.27 fps, 
p 0 = 0.07513 lb/ft 8 , 

M 0 = v 0 /a, (64) 

M = v cos -S/a cos $ ot (65 ) 

Y = poid 2 /m. (66) 

Then the horizontal range in feet is 

x =(Po/PY)[S(M) -S(M„)], (67) 

the time of flight in seconds is 

t f = (p 0 sec VapY) [T(M) - T(M.) ], (68) 
the altitude in feet is 


y = x tan # 0 


tan S = tan S 0 


aop 0 2 sec 2 &o 

a 2 p 2 Y 2 

^Po sec2 -S 0 jj(UJ) 


a 2 PY 


(A(M)-A(M.)] +x 
I(M 0 ) ]. 


appo sec 2 $o 
a 2 PY 


I(M 0 ). 


The angle of departure for y = 0 may be found by the formula 

sin 2S a = . A (^) ~7 A ( M °) . _ j(M 0 ) j . 

a 2 PYL SfM) -SCMJ J 


(69) 

(70) 

(71) 


If p = po and a = a<>, let S = S(M), S 0 = 
S(M 0 ), etc. Then these formulas can be ex- 
pressed more simply: 

x = (1 /y) (S — S 0 ), (67a) 

t f = (sec. •So/aoY) (T — T,), (68a) 


y = x tan S 0 


sec 2 S 0 
a* y 2 


(A — A.) + x 


sec 2 S 0 T 
aoY Io 
(69a) 


tan S = tan S 0 
and for y = 0, 

2 

sin 2S 0 = 

a«Y 


sec 2 S 0 

a* Y 


(I-Io), 


r A -Ao i i 
Ls-So lo J 


(70a) 

(71a) 
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(b) Example 


Data: Gun, 90-mm; Shell, HEAT, T108E20 


Caliber 

Mass 

Muzzle Velocity 
Form Factor 
Air Density 
Sound Velocity 


0.29S25 ft 
14.2 lb 
2400 fps 
1.00 on Gtios 
0.07513 lb/ft 3 
1120.27 fps 


Problem: To find the velocity, time of flight, angle 
of departure, and angle of fall at a range 
of 3000 feet. 


Solution : From the Table of the Siacci Functions 
based on Shell, 90-mm, HEAT, T108. 


Check: x tan 3 0 = 3000 (.00964) =28.92 ft 
sec 2 £ 0 (A -A 0 ) _ 1.00010(0.01867) 

aoT 2 ~ 0.000,238.29 

= -78.36 ft 

*^ 30000 ^ 0 . 00851 ) 

aoif 0.51667 

= 4. 49.42 ft 

y = 28.92 + 49.42 - 78.36 = — 0.02 ft 
(Eq. 69a). 

7. Differential Effects 
Some of the effects on range and deflection 
caused by variations from standard conditions 
will be discussed and explained. More detailed 
explanations are given in References 21, 22, 23 
and 24. 


MS T A I 

1.73650 3.091% 1.427502 .02537 .01923 

2.14234* 1.70836 0.709983 .00670 .00851 

Dif! 1.38360** 0.717519 .01867 .01072 

* M. = 2400/1120.27 = 2.14234 (Eq. 64) 

T = 0.07513(1.00) (0.29525 ) 714.2 = 
0.000,461,2 ft- 1 (Eq. 66) 

aoY = 1120.27(0.000,461,2) = 0.51667 sec" 1 

ao,Y a = 0.51667(0.000,461,2) = 

0.000,238.29 ft" 1 sec" 1 

**S — So = -pc = 0.000,461,2(3000) = 

1.38360 (Eq. 67a) 

si " 24 -= 03 w[tS- 0 00851 ] 

= 0.01928 (Eq. 71a) 

2# 0 .= 19.64 mils, £ 0 = 9.82 mils, cos $0 = 
0.99995, 

sec = 1.00005, sec 2 $ 0 = 1.00010 

tan = 0.00964 

1 OOOIO 

- 0 51567 (0*01072) = - 0.02075 

tan £ = -0.01111 (Eq. 70a) 

$ = — 1 1.32 mils, cos & = 0.9994 
v = 1120.27 (1.73650) (0.99995)/(0.99994) 

= 1945.37 fps (Eq. 65) 

1 OOOO s 

t, = (0.717519) = 1.389 sec (Eq. 68a) 


a. Effects on Range 

(1) Height 0 } Target . Two methods may be 
used for determining the effect of the height of 
the target with respect to the gun : one by means 
of the trajectory chart, and another by a correction 
to the angle of departure. 

(a) For large variations in height, the tra- 
jectory chart should be used. The trajectory that 
passes through the point determined by the given 
horizontal range and altitude can be found. The 
total range of this trajectory is the range where 
it crosses the axis. The difference between the 
given range and the total range is the effect of the 
height of target. For greater accuracy, the ex- 
terior ballistic tables may be used in lieu of the 
charts. 

(b) For small variations in height, an ap- 
proximate effect may be found by adding the 
angle of site, 

« = tan“ 1 (H/Ro), (72) 

to the angle of departure 9 0 , obtaining the cor- 
rected angle 

<P = ?o + e. (73) 

Here, H is the height of target, 

R<> — the map range, 

9o — the angle of departure corresponding to 
the map range in the Exterior Ballistic 
Tables, 

? — the angle of departure corresponding to the 
corrected range R. 

Then the effect on range due to height of target is 
A h R — R<> — R. (74) 
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(b) Example 


Data: Gun, 90-mm; Shell, HEAT, T108E20 


Caliber 

Mass 

Muzzle Velocity 
Form Factor 
Air Density 
Sound Velocity 


0.29S25 ft 
14.2 lb 
2400 fps 
1.00 on Gtios 
0.07513 lb/ft 3 
1120.27 fps 


Problem: To find the velocity, time of flight, angle 
of departure, and angle of fall at a range 
of 3000 feet. 


Solution : From the Table of the Siacci Functions 
based on Shell, 90-mm, HEAT, T108. 


Check: x tan 3 0 = 3000 (.00964) =28.92 ft 
sec 2 £ 0 (A -A 0 ) _ 1.00010(0.01867) 

aoT 2 ~ 0.000,238.29 

= -78.36 ft 

aoif 0.51667 

= 4. 49.42 ft 

y = 28.92 + 49.42 - 78.36 = — 0.02 ft 
(Eq. 69a). 

7. Differential Effects 
Some of the effects on range and deflection 
caused by variations from standard conditions 
will be discussed and explained. More detailed 
explanations are given in References 21, 22, 23 
and 24. 


MS T A I 

1.73650 3.091% 1.427502 .02537 .01923 

2.14234* 1.70836 0.709983 .00670 .00851 

Dif! 1.38360** 0.717519 .01867 .01072 

* M. = 2400/1120.27 = 2.14234 (Eq. 64) 

T = 0.07513(1.00) (0.29525 ) 714.2 = 
0.000,461,2 ft- 1 (Eq. 66) 

a.? = 1120.27(0.000,461,2) = 0.51667 sec" 1 

ao,Y a = 0.51667(0.000,461,2) = 

0.000,238.29 ft" 1 sec" 1 

**S — So = -pc = 0.000,461,2(3000) = 

1.38360 (Eq. 67a) 

si " 24 -= 03 w[tS- 0 00851 ] 

= 0.01928 (Eq. 71a) 

2# 0 = 19.64 mils, £ 0 — 9.82 mils, cos $0 = 
0.99995, 

sec = 1.00005, sec 2 $ 0 = 1.00010 

tan = 0.00964 

1 OOOIO 

- 0 51567 (0*01072) = - 0,02075 

tan £ = -0.01111 (Eq. 70a) 

$ = — 1 1.32 mils, cos & = 0.9994 
v = 1120.27 (1.73650) (0.99995)/(0.99994) 

= 1945.37 fps (Eq. 65) 

1 OOOO s 

t, = (0.717519) = 1.389 sec (Eq. 68a) 


a. Effects on Range 

(1) Height 0 } Target . Two methods may be 
used for determining the effect of the height of 
the target with respect to the gun : one by means 
of the trajectory chart, and another by a correction 
to the angle of departure. 

(a) For large variations in height, the tra- 
jectory chart should be used. The trajectory that 
passes through the point determined by the given 
horizontal range and altitude can be found. The 
total range of this trajectory is the range where 
it crosses the axis. The difference between the 
given range and the total range is the effect of the 
height of target. For greater accuracy, the ex- 
terior ballistic tables may be used in lieu of the 
charts. 

(b) For small variations in height, an ap- 
proximate effect may be found by adding the 
angle of site, 

« = tan“ 1 (H/Ro), (72) 

to the angle of departure 9 0 , obtaining the cor- 
rected angle 

<P = ?o + e. (73) 

Here, H is the height of target, 

R<> — the map range, 

9o — the angle of departure corresponding to 
the map range in the Exterior Ballistic 
Tables, 

? — the angle of departure corresponding to the 
corrected range R. 

Then the effect on range due to height of target is 
A h R — R<> — R. (74) 
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(2) Elevation . Since the Exterior Ballistic 
Tables give the range as a function of elevation for 
constant values of muzzle velocity and ballistic 
coefficient, the effect of an increase in elevation 
may be found from the differences. At short 
ranges, the effect may be computed more ac- 
curately by using the Siacci tables, as explained 
above. The vacuum formula (Eq. 6) may be use- 
ful in estimating the effect, especially for large 
ballistic coefficients. It should be noted that the 
effect is small in the vicinity of the maximum 
range and, of course, is negative at elevations 
above that corresponding to the maximum range. 

(3) Muzzle Velocity. Since the Exterior Bal- 
listic Tables give the range as a function of muzzle 
velocity for constant values of elevation and 
ballistic coefficient, the effect of an increase in 
muzzle velocity may be found from the differences. 
At short ranges, the effect may be computed by 
using the Siacci tables. In vacuum, the range is 
proportional to the square of the velocity, as shown 
in Equation 6 ; in air, this is a good approximation 
if the ballistic coefficient is large. 

(4) Ballistic Coefficient . Since the Exterior 
Ballistic Tables give the range as a function of 
ballistic coefficient for constant values of muzzle 
velocity and elevation, the effect of an increase in 
log C may be found from the differences. At 
short ranges, the effect may be computed by using 
the Siacci tables. It is customary to find the effect 
of a 1% increase in C, which is equivalent to an 
increase of 0.00432 in log C. 

(5) Weight of Projectile . An increase in 
weight of projectile decreases the muzzle velocity 
and increases the ballistic coefficient. 

The effect on muzzle velocity may be expressed 
ApV 0 = nv 0 (Ap)/p, (7$) 

where v 0 is the muzzle velocity 
p is the projectile weight, 
n the logarithmic differential coefficient. 

The value of n may be obtained from interior 
ballistic tables. It is approximately — 0.3 for a 
rifled gun with multiperforated propellant grains, 
— 0.4 for a rifled gun with single-perforated 
grains, — 0.47 for a smooth bore mortar with 
flake propellant, and — 0.6S for a recoilless rifle 
with multiperforated grains. 

Since the ballistic coefficient is proportional to 
the projectile weight by definition. Equation (22), 


the effect on C may be expressed 

A p C = C(Ap)/p. (76) 

If A V R is the effect on range of 1 unit increase in 
muzzle velocity, AcR the effect on range of 1% 
increase in ballistic coefficient, the effect on range 
of an increase in projectile weight is 

A P R = A v R(Apv) + A C R(100 A P C)/C. (77) 

(6) Air Density. The retardation due to air 
resistance is proportional to the air density and 
inversely proportional to the ballistic coefficient, 
as indicated in Eq. (25). Therefore, to a first 
approximation, the effect on range of an increase 
in air density is equal to the effect of the same 
percentage decrease in C. 

(7) Air Temperature . An increase in the 
temperature of the atmosphere affects the range in 
two ways: it increases the density of the air and 
the velocity of sound. The temperature of the 
air may affect the temperature of the propellant, 
which affects the muzzle velocity, but that is an 
interior ballistic problem. 

For a given pressure and humidity, the density 
is proportional to the absolute temperature. The 
effect of density on range has been treated above. 

For a given humidity, the velocity of sound is 
proportional to the square root of the absolute 
temperature. For a given projectile velocity, the 
Mach number is inversely proportional to the 
acoustic velocity. A change in M produces a 
change in the drag function. A new trajectory 
could be computed with the varied function ; how- 
ever, it is simpler to use an approximate formula 
deduced from a semi-empirical formula derived by 
Dr. Gronwall: 36 

A t R = 0.1929 [0.01R - 0.005v o A v R 
- (1 - 0.000, 015y s )A o R] At, (78) 

where R is the range, 

A t R~ the increase in range due to increase of 
temperature 

v 0 — the muzzle velocity in feet per second, 

A V R — the effect on R of 1 fps increase in v 0 , 

y B — the maximum ordinate in feet, 

A c R — the effect on R of 1% increase in C, 

At — the increase in temperature in degrees 
Fahrenheit. 
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( 8 ) Wind . A rear wind has three effects: it 
moves the coordinate system along with it rela- 
tive to the earth; it increases the angle of de- 
parture relative to the air; and it decreases the 
muzzle velocity relative to the air. The total 
effect on range may be expressed as follows : 

ao ur Tt , si n 9 A 9 R A V R1 

w =wi |_ T +— -ar- ros? 'sd- 

(79) 

where W R is the range component of the wind (a 
rear wind is positive; a head wind is negative), 

T — the time of flight, 

9 — the angle of departure, 

v 0 — the muzzle velocity, 

A 9 R— the effect on range due to a change of 
A? radians in 9 , 

A*R — the effect on range due to a change Av 0 in 

V 0 . 

b . Effects on Deflection 

(1) Wind. A crosswind has two effects: it 
moves the coordinate system with it and it changes 
the direction of motion relative to the air. The 
total effect on linear deflection may be expressed 

A W Z = W, (T — R/v 0 cos 9 ), (80) 

where W, is the crosswind. 

The angular deflection due to wind, in radians, 
is 


D w = W, (T/R — l/v 0 cos 9 )*. (81) 


(2) Drift. After the yaw due to initial dis- 

turbances has damped out, a spinning shell has 
a yaw of repose due to the curvature of {he tra- 
jectory. This yaw is mainly to the right of the 

trajectory if the shell has right-handed spin; to 
the left if it has left-handed spin. It* may be de- 

termined at short ranges by firing alternately 

from barrels with right and left hand twists of 

rifling, or at long ranges by correcting the ob- 

served deflection for the effects of wind, cant of 
trunnions, and rotation of earth. 


The differential equation of motion for the drift 
z is 28 


.. _2xgAvoLN . Gpz 
z mdnu*MN 0 X p 0 C ’ 


(82) 


where g is the gravitational acceleration, 

A — the axial moment of inertia, 
v 0 — the muzzle velocity, 

L — the crosswind force (lift), 

N — the spin, 

N 0 — the initial spin, 
m — the mass of the projectile, 
d — the caliber, 

u — the velocity of the projectile relative to the 
air, 

M — the overturning moment about the center of 
gravity, 

x — the horizontal range, 

G — the drag function, 
p — the air density, 

p 0 — the standard air density (0.07513 lb/ft 2 * * * * * 8 ), 

C — the ballistic coefficient. 


Dots denote time derivatives, and the initial con- 
ditions are 

Zo = 0, Zo = 0. 

c. Ballistic Density , Temperature, and Wind 

In treating the effects of variations in the 
atmosphere, it has been assumed that the percent- 
age variations from standard density and tempera- 
ture and the wind are constant all along the tra- 
jectory. Actually, they vary with altitude, so 
that weighted mean values must be used. The 
weighting factors are determined by the propor- 
tion of the total effect that a variation occurring 
in each altitude zone would have. Instead of 
computing a set of weighting factors for each tra- 
jectory, a few representative curves are used. The 
sum of the products of the zonal variations by 
their weighting factors is called the ballistic 
density, temperature, or wind. The results are 
then treated as constants in computing the effects 
on range and deflection. 
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GLOSSARY 


angle of departure, (elevation) As used in this 
pamphlet these terms mean quadrant angle of 
departure which is defined as the angle between 
the horizontal and the line along which the pro- 
jectile leaves the projector. One of the two 
initial conditions of a trajectory, the other con- 
dition being initial velocity . 

For a rocket-fired projectile in which the pro- 
pellant burning continues after the rocket leaves 
the launcher, it is necessary to make assump- 
tions relative to the point at which the rocket 
becomes a missile in free flight and the velocity 
attained at that point. 

angle of inclination. The angle between the hori- 
zontal and the tangent to the trajectory at any 
point. At the origin of the trajectory this is the 
same as the angle of departure . 

apical (semi-apical) angle. In general the angle 
formed at the apex or tip of anything. As ap- 
plied to projectiles, the angle between the tan- 
gents to the curve outlining the contour of the 
projectile at its tip, or for semi-apical angle, 
the angle between the axis and one of the tan- 
gents. For a projectile having a conical tip, 
the cone apex angle. 

ballistics. Branch of applied mechanics which deals 
with the motion and behavior characteristics of 
missiles, that is, projectiles, bombs, rockets, 
guided missiles, etc., and of accompanying phe- 
nomena. It can be conveniently divided into 
three branches: — interior ballistics, which deals 
with the motion of the projectile in the bore of 
the weapon ; exterior ballistics, which deals with 
the motion of the projectile while in flight; and 
terminal ballistics, which is concerned with the 
effect and action of the projectile when it im- 
pacts or bursts. 

crosswind force. The component of the force of 
air resistance acting in a direction perpendicular 
to the direction of motion, and in the plane de- 
termined by the tangent to the trajectory formed 
by the center of gravity of a missile and the axis 
of the missile. (This plane is called the plane 
of yaw.) 

drag. The component of the force of air resistance 
acting in a direction opposite to the motion of 
the center of gravity of a missile. 


elevation. See angle of departure. 

exterior ballistics. See ballistics . 

horizontal range. The horizontal component of 
the distance from the beginning of the trajectory 
(point of origin) to any other point on the tra- 
jectory. 

initial velocity, (muzzle velocity) The projectile 
velocity at the moment that the projectile ceases 
to be acted upon by propelling forces. 

For a gun-fired projectile the initial velocity, 
expressed in feet or meters, per second, is called 
muzzle velocity or initial velocity . It is obtained 
by measuring the velocity over a distance for- 
ward of the gun, and correcting back to the 
muzzle for the retardation in flight. 

For a rocket-fired projectile, the velocity at- 
tained at the moment when rocket burnout oc- 
curs is the initial velocity. A slightly fictitious 
value is used, referred to an assumed point as 
origin of the trajectory. 

The initial velocity of a bomb dropped from 
an airplane is the speed of the airplane. 

Magnus force. The component of the force of air 
resistance acting in a direction perpendicular to 
the plane of yaw (see crosswind force). It is 
caused by the action between the boundary layer 
of air rotating with the projectile and the air 
stream. This force is small, compared with 
drag and crosswind force . 

muzzle velocity. The velocity with which the 
projectile leaves the muzzle of the gun. See 
initial velocity . 

slant range. The distance, measured along a 
straight line, from the beginning of the trajectory 
(point of origin) to a point on the trajectory. 

stability factor. A factor which indicates the rela- 
tive stability of a projectile under given condi- 
tions. It is dependent upon the force moments 
acting to align the projectile axis with the tan- 
gent to the trajectory, and the overturning force 
moments. For the projectile to be stable, the 
stability factor must be greater than unity. 

standard atmosphere. An average or representa- 
tive structure of the air used as a point of de- 
parture in computing firing tables for projec- 
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tiles. The standard atmosphere adopted in 1956 
for the Armed Services and now known as the 
U. S. Standard Atmosphere, is that used by the 
International Civil Aviation Organization 
(ICAO). This Standard Atmosphere assumes 
a ground pressure of 760 m/m of mercury, a 
ground temperature of 15° C, air density at 
ground level of 0.076,475 lb/ft 3 , and a tempera- 
ture gradient to the beginning of the strato- 
sphere expressed by the formula 

absolute temperature T( C K) = 288.16-6.5H 

where H is height above sea level in kilometers. 

In the stratosphere (above 11 kilometers) the 
temperature is assumed constant at 216.66° K. 
For firing table purposes wind velocities of zero 
are assumed. 


superelevation. An added positive angle in anti- 
aircraft gunnery that compensates for the fall of 
the projectile during the time of flight due to the 
pull of gravity; the angle the gun or projector 
must be elevated above the gun -target line. 

trajectory. The path followed in space by the 
center of gravity of a missile while in free flight. 

yaw. The angle between the direction of motion 
of a projectile and the axis of the projectile, re- 
ferred to either as yaw , or more completely, 
angle of yaw . The angle of yaw increases with 
time of flight in an unstable projectile (see sta- 
bility factor ) and decreases to a constant value, 
called the yaw of repose, or the repose angle of 
yaw, in a stable projectile. 
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Figure 2. Drag Coefficient vs Mach Number 
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Figure 3. Time of Flight for Projectile Type 1 
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Figure 6 (1 of 12). Trajectories for Projectile Type 2, 45° elevation, 1200 fps muzzle velocity 




Figure 6 (2 of 12). Trajectories for Projectile Type 2, 45° elevation, 2000 fps muzzle velocity 












Figure 6 (3 of 12). Trajectories for Projectile Type 2, 45° elevation, 2800 fps muzzle velocity 
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Figure 6 (5 of 12). Trajectories for Projectile Type 2, 45° elevation, 4400 fps muzzle velocity 
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Figure 6 (6 of 12). Trajectories for Projectile Type 2, 45° elevation, 5200 fps muzzle velocity 





Figure 6 (7 of 12), Trajectories for Projectile Type 2, 60° elevation, 1200 fps muzzle velocity 
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Figure 6 (8 of 12). Trajectories for Projectile Type 2, 60° elevation, 2000 fps muzzle velocity 
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Figure 6 (9 of .12). Trajectories for Projectile Type 2, 60® elevation, 2800 fps muzzle velocity 
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Figure 6 (12 of 12). Trajectories for Projectile Type 2, 60° elevation, 5200 fps muzzle velocity 












